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Abstract. This chapter takes a microscopic view of quantum tunneling of mag- 
netization (QTM) in single-molcculc magnets (SMMs), focusing on the interplay 
between exchange and anisotropy. Careful consideration is given to the relation- 
ship between molecular symmetry and the symmetry of the spin Hamiltonian that 
dictates QTM selection rules. Higher order interactions that can modify the usual 
selection rules are shown to be very sensitive to the exchange strength. In the strong 
coupling limit, the spin Hamiltonian possess rigorous D2h symmetry (or Coo in high- 
symmetry cases). In the case of weaker exchange, additional symmetries may emerge 
through mixing of excited spin states into the ground state. Group theoretic argu- 
ments arc introduced to support these ideas, as arc cxtcuisivc results of magnetization 
hysteresis and electron paramagnetic resonance measurements. 



1.1 Spin Hamiltonian 

The concept of an effective spin-Hamiltonian involving only spin variables 
has been employed in the study of paramagnetic species for well over half a 
century. This formalism is particularly suited to the study of transition metal 
complexes in vi^hich the groimcl state is very often an orbital singlet that is 
well isolated from excited orbital states due to the strong influence of the 
ligand field [1]. The ground state multiplicity is determined entirely by the 
spin state of the ion in this situation, even though the ground wave functions 
are not exact eigenstates of because of residual spin-orbit (SO) coupling. 
It is this coupling that gives rise to the familiar anisotropic zero-field splitting 
(zfs) and Zeeman terms in the resultant spin Hamiltonian, parameterized by 

the zfs D and Lande ^-tensors. 
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1.1.1 Giant-Spin Approximation Hamiltonian 

The magnetic moment of a typical polynuclear transition metal cluster is de- 
termined by the exchange interactions between the spins associated with the 
constituent ions. As detailed in this chapter, there arc a number of ways to 
extend the spin Hamiltonian formalism to this multi-ion situation. By far 
the simplest is the so-called Giant Spin Approximation (GSA), in which one 
assigns a total (giant) spin quantum number, S, to the lowest-lying (m^) 
magnetic levels [2]; for a ferromagnetic molecule, S is obtained from the alge- 
braic sum of the spin values associated with each of the ions. If the exchange 
coupling within the molecule is large in comparison to the single-ion zfs inter- 
actions, then this ground spin multiplet will be well separated from excited 
spin states. One may then employ a GSA Hamiltonian to describe the mag- 
netic properties of the molecule, provided that the temperature is sufficiently 
low that excited spin states are not thermally populated. 

A series expansion in terms of the spin component operators S^, Sy, and 
Sz, employing so-called Extended Stevens operators, results in the following 
effective zfs Hamiltonian [3-5]: 

2S p 

^^f^ = EE^p^P' (1-1) 

p q=0 

where 0'^{Sx, Sy, S^) represent the operators, and the associated phc- 
nomenological (or effective) zfs parameters. The subscript, p, denotes the order 
of the operator, which must be even due to the time reversal invariance of the 
SO interaction; the order is also limited by the total spin, S, of the molecule 
such that p < 2S. The superscript, q {< p), denotes the rotational symmetry 
of the operator about the ^-axis. Eqn. 1.1 has been employed with great suc- 
cess in the study of single-molecule magnets (SMMs), particularly in terms of 
describing low-temperature quantum tunneling of magnetization (QTM) be- 
havior and electron paramagnetic resonance (EPR) data [2]. In fact, Eqn. 1.1 
has even been applied quite successfully in cases where the ground spin mul- 
tiplet is not so well isolated from excited spin states [6-9]. Fourth and higher 
order operators are often found to be important in these cases. This chapter 
examines the microscopic origin of these terms, which are often negligible for 
the constituent ions. However, the usually dominant 2"^-order zfs interaction 
is first considered. 

Second Order Anisotropy 

Although fourth and higher order Stevens operators (p > 4) are allowed for a 
molecule with 5 > 2, it is not always necessary to include all of them (up to 
p = 25). For SMMs comprised of transition metal ions, the low-energy physics 
is usually dominated by 2"'^-order SO anisotropy. Therefore, the simplest zero- 
field GSA Hamiltonian used in characterizing SMMs is often written as: 
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H,,, = DSl + E{Sl-Sl). (1.2) 

Eqn. 1.2 includes only 2"'*-order terms, where D (= S-B^) parameterizes the 
uniaxial anisotropy and E (= B^) the rhombicity. For an approximately uniax- 
ial system, DS^ is the dominant anisotropy, with z chosen as the quantization 
axis. In biaxial cases, the ratio between E and D is usually restricted such 
that \E/D\ < 1/3; one can always perform a rotation of the coordinate system 
such that this criterion is satisfied. 

One of the main goals of this section is to understand the influence of 
molecular symmetry on the QTM properties of SMMs. Hence, it is important 
to examine the symmetry of Eqn. 1.2 since, strictly speaking, the symmetry 
of the Hamiltonian should be compatible with the symmetry of the molecule 
under investigation. In addition, the nature of SO coupling ensures that the 
spin Hamiltonian be invariant under time-reversal {p is even), i.e., the spin 
Hamiltonian naturally possesses Cj symmetry. As such, the physics is invariant 
to inversion of either the total spin moment or the applied field. A classical 
representation of Eqn. 1.2 is shown in Fig. 1.1, with \E\ = \D\/5 and D < 0. 
This graphical representation is obtained by substituting the spin operators 
in Eqn. 1.2 by their classical equivalents, as follows: 

Sx ^ S sin 6 cos (f); 
Sy ^ S sin 9 sin 0; 
Sz ^ S cos 6; 

where and (f> arc the inclination and azimiithal angles in spherical coordi- 
nates, respectively. In this representation, the spin is treated as a macroscopic 
magnetic moment for which all the three components (Sx, Sy and Sz) can 
be determined simultaneously. The surface shown in Fig. 1.1 represents the 
energy of the spin as a function of its orientation, where the radial distance 
to the surface corresponds to its energy. As can be seen, Eqn. 1.2 contains the 
following symmetry elements: (a) three orthogonal C2 axes, corresponding to 
X, y and z, and (b) three orthogonal mirror planes, corresponding to the xy, 
yz, and 2a;-planes. These symmetry elements, together with the Cj symmetry, 
give rise to a -D2/1 symmetry for Eqn. 1.2, which is obviously a much higher 
symmetry than most real molecules. Consequently, even though Eqn. 1.2 can 
often account very well for low-temperature thermodynamic measurements 
performed on SMMs {e.g. ac susceptibility and magnetization), it may never- 
theless fail to explain symmetry-sensitive quantum mechanical (spectroscopic) 
observables such as QTM steps. Berry-phase interference (BPI) patterns and 
EPR/neutron spectra. 

It should be noted from the preceding discussion that the B2O2 {O^ = 
^[SxSz + SzSx]) term was neglected in Eqn. 1.2, even though it is perfectly 
allowed within the GSA. With the exception of Cj, the 0| operator satisfies 



(1.3) 
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none of the symmetry operations described in the preceding paragraph. How- 
ever, such a term is unnecessary, as can be seen when writing the 2"'^-order 
GSA Hamiltonian in the more compact form: 

= S-D S, (1.4) 

where Z) is a 3 x 3 matrix corresponding to the full 2"''-order anisotropy tensor. 

In Eqn. 1.4, D and E are related to the diagonal elements of D (see below) 

while B2 appears as off-diagonal elements. The only restriction on D is that 
it must be Hermitian in order to guarantee the Hamiltonian be Hcrmitian; 

indeed, D^z = D^x = ^B^. Consequently, D can always be diagonalized by 
rotating the original Cartesian coordinate frame. Upon doing so, all of the 
off-diagonal elements of the rotated matrix vanish, i.e.. B\ = in the new 
Cartesian coordinate frame. Finally, one may adjust the absolute values of 
the resultant eigenvalues without altering the symmetry of the Hamiltonian 

simply by subtracting ^{D^x + Dyy)I from _D (/ is the identity matrix). The 
zfs Hamiltonian can then be rewritten as Eqn. 1.2 with 

D = - ^{Dxx + Dyy) and E = ^{Dxx - Dyy), (1.5) 

where Da (i = x, y, z) refer to components of the diagonalized (rotated) D 
tensor. In other words, Eqn. 1.4 is equivalent to Eqn. 1.2, requiring just two 
parameters, D and E, to completely describe the effective 2"'^-order anisotropy 
within the GSA. Inclusion of 0\ results simply in a rotation of the surface 
depicted in Fig. 1.1. Consequently, the 2"^-order GSA Hamiltonian necessarily 
possesses at least symmetry. 

The preceding discussion demonstrates something very important: even 
though 2'^'^-order terms typically represent the dominant interactions within 

the GSA description of a SMM, the resultant D tensor possesses an artifi- 
cially high {D2h) symmetry which may not be compatible with the structural 
symmetry of a particular molecule under investigation. The consequences of 
this property of the GSA in terms of the resultant QTM will be discussed in 
detail in following sections. 

Higher-Order Anisotropies 

For a SMM with S > 2, Stevens operators of order four (and higher) are 
allowed in the GSA Hamiltonian— note p can take on any even value from 
2 to 25". The values of the 4*'^-ordcr parameters arc often deceptively small, 
especially for SMMs with large spin values. For example, \B^/D\ ~ 5 x 10~^ 
for the Mni2 SMMs, yet the B^O^ GSA term contributes ~ 20% to the energy 
barrier. This is due not only to the higher order of Sz in O4, but also because 
of the way in which the O4 operator is defined— a multiplier of 35 is associated 
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with S'^. In general, the contribution of higher-order terms to the energies of 
spin states may be expected to be smaller than those of the 2"'^-ordcr terms. 
However, this rule of thumb breaks down in the weak exchange limit (or for 
particularly high-symmetry molecules [10]); indeed, it is in this limit that 
one may call into question the validity and/or usefulness of the GSA. Axial 
(g = 0) 4*^-order terms lead to a non-parabolic energy barrier, which gives rise 
to non-even spacings between EPR and QTM resonance fields [11-14]. More 
importantly, the higher-order transverse {q ^ 0) terms introduce additional 
symmetries into the GSA Hamiltonian, enabling a more precise description of 
the quantum properties of SMMs. 

Fig. 1.2 displays the classical energy surfaces corresponding to the 4**^- 
order Stevens operators; the 0% surface is not shown since it commutes with 
Sz and possesses Coo (cylindrical) rotational symmetry. All of the surfaces, 
and hence the operators, exhibit rotational symmetries which are compatible 
with the superscript q. However, one may note a systematic difference between 
the g-odd and g-even operators: the g-even operators have the a;y-plane as an 
additional mirror plane; this is not the case for the g-odd operators. In par- 
ticular, the lobes of the 0| and 0\ operators lie within the a;2/-plane, whereas 
those of 0\ and 0\ alternately lie above and below this plane. The symme- 
tries of the 4"^-ordcr operators can be understood in terms of combinations of 
rotational symmetries and the intrinsic Cj symmetry of the Hamiltonian. For 
the g-even operators, the direct product of the rotational and inversion group 
leads to: C2 x Ci = C2/1 symmetry for 0|; and C4 x = Ca^u for 0|. Thus, the 
xy-plane is introduced as a new symmetry element. In contrast, for the g-odd 
operators, C\ x Ci = Ci for O4 and C3 x Ci = Sq for Of. The resultant sym- 
metry groups corresponding to these operators include an improper rotation 
{C'i can be treated as the improper rotation 82)- The absence of the xy-mirror 
plane for the g-odd operators suggests that the molecular hard plane may 
not coincide with the xy-plane, which leads to several intriguing phenomena 
described later in this chapter. 

The inclusion of p > 4 Stevens operators in the GSA has a significant influ- 
ence on the interpretation of QTM measurements. When limited to 2"'^-order 
anisotropy, the zero-field Hamiltonian can only mix spin projection states that 
differ in by an even number, i.e., ArUg = |msi — ms2\ = 2n (n =integer). 
This means that only fc-cven (k = rrigi +771^2) QTM resonances should be ob- 
servable for parallel applied fields {H || z). Moreover, in molecules for which 
rhombicity is symmetry forbidden (£^ = 0), a purely 2"'^-order Hamiltonian 
would be cylindrically symmetric. Consequently, for H \\ z, remains an 
exact quantum number and QTM should be completely forbidden. However, 
plenty of such high-symmetry molecules exist and are known to exhibit clear 
QTM behavior, including Mni2 and several other SMMs discussed in this 
chapter. In these situations, it is necessary to include higher order anisotropies 
in the GSA. The corresponding operators introduce new spin-mixing rules 
which can lead, for example, to fc-odd QTM resonances. 
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The advantage of the GSA lies in the fact that one can usually restrict 
the total number of zfs parameters involved in data analysis to just a few by 
considering the overall symmetry of the molecule under study. Furthermore, 
the GSA Hilbert space includes only the 25+1 states that belong to the 
ground spin multiplct, such that the Hamiltonian matrix has dimension {2S + 
1) X {2S + 1). This makes data analysis for large clusters computationally 
possible. However, the GSA completely ignores the internal degrees of freedom 
within a molecule, thus completely failing to capture the underlying physics 
in cases where the total spin fluctuates [15-17]. Moreover, when a molecule 
possesses very little symmetry {e.g. Ci), the number of GSA zfs parameters 
cannot be restricted on the basis of symmetry and, in principle, all possible 
terms (up to p = 2S) should be taken into account. In these cases, it may be 
advantageous to employ a multi-spin Hamiltonian, particularly in situations 
where microscopic insights are desired. 

1.1.2 Multi-Spin Hamiltonian 

In the multi-spin (MS) model, a molecule is treated as a cluster of magnetic 
ions (spins) which are coupled to each other via pairwise exchange interactions. 
The corresponding zero-field Hamiltonian is: 



where Sj represents the spin operator of the i^^ ion, and di is the 2"'^-order 
zfs tensor associated with this same ion: lowercase symbols arc used here to 
differentiate parameters/variables employed in both models, i.e., lowercase = 

MS and uppercase = GSA. For the sake of simplicity, di is written in the 
diagonal form: 



where the local coordinate frame of di is chosen to match the local principal 
anisotropy axes of the i ion. Ri is the Euler matrix, specified by the Euler 
angles 9i, (pi and ipi, which transforms the local coordinate frame of the i"^ ion 

into the molecular coordinate frame. The matrix Jij specifies the exchange 
interaction between the i^^ and j**^ ions. It should be emphasized that all 
of the parameters in Eqn. 1.6 should be constrained by the structure of the 
molecule under study, i.e., the overall symmetry of the Hamiltonian must be 
compatible with the molecular symmetry. 

The MS model captures physics associated with internal molecular degrees 
of freedom that are not easily understood within the GSA framework. First 




(1.6) 



i i<j 



Cj 
d = -e^ 




(1.7) 
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and foremost, the MS model is capable of describing phenomena in which the 
total spin of a molecule fluctuates, i.e., it gives the energies of excited spin 
states in addition to the ground state, and includes the mixing between these 
states [18]. Secondly, the parameters in the MS Hamiltonian have clear phys- 
ical significance, i.e., they describe the magnetic properties of the constituent 
ions and the coupling between them. Moreover, many of these parameters can 
be independently verified through measurements of related compounds [19]. 
In contrast, the parameters deduced on the basis of a GSA arc purely phe- 
nomenological. For example, comparisons between the two models have shown 
that higher order anisotropics in the GSA arise from the interplay between the 
local 2"'^-order single-ion anisotropy and the magnetic interactions between 
the ions, leading to mixing of excited spin states into the ground spin multi- 
plet [16,18,20]. In other words, the phenomenological p > 4 zfs parameters 
are a direct manifestation of physics that goes beyond the GSA. The following 
section deals with this issue in detail, with a focus on the correlation between 
QTM behavior and the structural symmetries of real molecules. 

1.2 Quantum Tunneling of Magnetization in 
High-Symmetry Mna Single-Molecule Magnets 

The first clear observation of QTM selection rules, i.e., a complete absence 
of symmetry forbidden resonances [21], was reported for the trigonal SMM 
[NE4]3[Mn3Zn2(salox)30(N3)6Cl2] (henceforth Mna) [22,23]. This section fo- 
cuses on QTM in SMMs with trigonal symmetry, emphasizing (i) symmetry- 
enforced selection rules that allow quantum relaxation in k-oAA resonances 
(Sec. 1.2.2), (ii) the role of disorder (Sec. 1.2.3), and (iii) the microscopic ori- 
gin of the B\0\ GSA interaction which is predicted to give rise to unusual 
BPI patterns (Sec. 1.2.4). 

1.2.1 The Mna Single-Molecule Magnet 

Several Mns SMMs are known to crystallize in the trigonal space group i?3c 
with racemic mixtures of C3 symmetric chiral molecules [18,22,23]. The struc- 
ture of the [NE4]3[Mn3Zn2(salox)30(N3)6Cl2] molecule is shown in Fig. 1.3a. 
The magnetic core consists of three ferromagnetically coupled Mn^^^ (s = 2) 
ions, which form an equilateral triangle with two Zn" ions located above and 
below the the Mn3 plane, thus forming a trigonal bipyramidal structure. The 
C3 axis of the molecule is perpendicular to the Mn3 plane, while the local 
easy-axes of the individual spins are defined by the Jahn- Teller (JT) elonga- 
tion axes of Mn^^^ ions, which are tilted slightly with respect to the C3 axis. 
At low temperatures, the spin 5 = 6 ground state multiplet can be described 
with the following GSA Hamiltonian: 



H = DSl + Biol + Biol + Biol + ,jlbB- g -S 



(1.8) 
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Due to the C3 symmetry of the molecule, the 2"'' order transverse anisotropy 
term, E^S"^ — Sy), is rigorously forbidden. Hence, the leading trigonal {0\) 

and hexagonal (Og) transverse zfs terms are instead included in Eqn. 1.8. 

Mn3 is highly attractive in the context of understanding the origin of QTM 
at a microscopic level. The dimension of the MS Hamiltonian matrix for three 
s = 2 spins is just [(2s + 1)^]^ = 125 x 125. The C3 symmetry reduces the 
number of interaction parameters to just a single exchange constant, J, and 
identical d and e values for each ion; it also guarantees identical 9i Euler an- 
gles (= 15°) for the three spins, with (f)i = — 1) x 120°. The remaining 
parameters have then been determined from fits to EPR and magnetization 
hysteresis measurements [18,21-23]. Lastly, the structure contains no solvent 
molecules. This is rare among SMMs and removes a major source of disor- 
der [24]. Consequently, exceptional spectroscopic data (QTM and EPR) are 
available against which one can test theoretical models. 

1.2.2 QTM Selection Rules in Mng 

For large spin systems, the effects oi q^Q zfs terms typically manifest them- 
selves at energy scales that are orders of magnitude smaller than those of the 
axial {q = 0) terms. One must therefore focus on avoided level crossings, where 
the tunneling gaps are governed by the transverse terms in Eqn. 1.8. Fig. 1.4 
displays the Zeeman diagram corresponding to the nominal spin 5 = 6 ground 
state multiplet of the Mns molecule. Due to symmetry restrictions {q = 3n 
for C3 symmetry, where n is an integer), non-zero tunneling gaps are limited 
to level crossings with Avris — 3n, where is the projection of the total spin 
onto the molecular C3 {z-) axis. All such gaps. A„i^.m/ , have been labeled in 
Fig. 1.4 for QTM resonances A; < 3, where k {— nis + m'^) denotes an avoided 
crossing between pairs of levels with spin projections and m'^ {fhs denotes 
-|m,|). 

By performing a mapping of the energy diagram obtained via exact diag- 
onalization of Eqn. 1.6 onto that of the GSA Hamiltonian (Eqn. 1.8) one can 
obtain microscopic insights into the emergence of p > 4 transverse terms in the 
latter approximation. Published zfs parameters were employed for simulations 
involving Eqn. 1.6, i.e., d = —4.2 K and e = 0.9 K [21]. An isotropic exchange 
constant J (= 10 K) was employed, set to a value that is artificially high in 
order to isolate the ground state from excited multiplets, thereby simplifying 
analysis of higher-lying QTM gaps (sec Fig. 1.4). The Euler angles were set 
to (pi ~ 0, (f)2 = 120° and (f>s = 240° (all t/ji = 0) to preserve C3 symmetry, 
while 0i (= 6) was allowed to vary in order to examine its influence on QTM 
selection rules. 

The situation in which the JT axes of the three Mn^^^ ions are parallel to 
the C3 axis is first considered, i.e., 6 = 0. The top section of Table. 1.1 lists 
the magnitudes of even-n QTM gaps involving pairs of levels with ArUg = 3n, 
deduced via diagonalization of Eqn. 1.6 in the absence of a transverse field, 
Ht (_L z). The odd-n, Ht = gaps are identically zero, as can be seen 
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from their dependence on Ht (Fig. 1.4 inset): the power-law behavior indi- 
cates no contribution from zfs interactions (at iJ-p = 0). Consequently, one 
expects only even-n zfs terms of the form i?p"Op" in the GSA: those satis- 
fying this requirement have six-fold rotational symmetry about the C3 axis, 
i.e., a higher symmetry than the real molecule (further explanation is given 
below). For comparison, these QTM gaps are simulated employing Eqn. 1.8 
with Bl=0 and = 4.3 x 10"^ K. As seen in Table 1.1, an excellent over- 
all agreement between the two models is obtained. Small differences may be 
attributed to higher-order six- fold terms such as BgOg, Bf^OfQ, etc., which 
have been neglected in this analysis. 

A more realistic situation involves a tilting of the JT axes away from 
the C3 axis hy = 8.5°, as is the case for Mna [23]. Both even- and odd-n, 
i?T = QTM gaps are generated in this situation, i.e., fc-odd QTM resonances 
become allowed. This may be understood within the framework of the GSA as 
being due to the emergence of zfs interactions possessing three-fold rotational 
symmetry about the molecular C3 axis, i.e., Bp^Op"^ with n = 1 and p > 
4; the leading such term is -B4O4. Table 1.1 lists the QTM gaps evaluated 
via diagonahzation of Eqn. 1.8 using i?| = 4.3 x 10"'^ K and = 4.77 x 
10~^ K. Excellent agreement is once again achieved between the GSA and MS 
Hamiltonians. Minor deviations may, in principle, be corrected by introducing 
higher-order transverse terms such as BqOq. 

The emergence of the BfOl interaction in the GSA description of Mn3 
clearly indicates a lowering of the symmetry of the spin Hamiltonian upon 
tilting the JT axes. To understand this one must consider both the symmetry 
of the molecule and the intrinsic symmetry of the zfs tensors of the individual 
ions. Considering only 2"*^ order SO anisotropy, the Hamiltonian of a single 
Mn''^ ion possesses symmetry (as noted in Sec. 1.1.1), with three mu- 
tually orthogonal C2 axes. When the JT axes are parallel {0 = 0), the local 
2-axis of each Mn^^^ center coincides with the molecular C3 axis. The resultant 
Hamiltonian should then possess C3 x C2 x Cj = Ceh symmetry (see Fig. 1.5a), 
requiring _B| = 0; the additional Cj symmetry arises from the time-reversal 
invariance of the SO interaction. In contrast, when the JT axes are tilted, the 
C2 and C3 axes do not coincide. In addition, the xy-mirror symmetry of the 
molecule is broken, as is that of the spin Hamiltonian. The rotational symme- 
try then reduces to three-fold and, hence, B4O4 is allowed; the symmetry in 
this case is C3 x Cj = Sq (Fig. 1.5b). 

It is possible to reinforce the preceding discussion via group theoretic ar- 
guments, without the need to write down an exact expression for the Hamil- 
tonian. When the external magnetic field is applied parallel to the molecular 
Cs-axis, the Ceh symmetry reduces to Cq, and the 13 basis functions of the 
5 = 6 ground state fall into six distinct one-dimensional irreducible repre- 
sentations [25]. These functions can be grouped according to their behavior 
under a Cg rotation: |-6) , |0) , \+6) € A; |-2) , \+4) G Ta; \+2) , |-4) e Tg; 
|-3) , \+3) e , |-5) e r^; j-l) ,|-|-5) e Te, where A.. .6 are the six irre- 
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ducible representations following the Bethe notation [25]. Because the Hamil- 

tonian operator belongs to the totally symmetric representation, {ms \ H \m'^) 
is non-zero only when \ms) and \m'g) belong to the same representation [26]. 
As can be seen, such states have Amg = 3n, with n even, which is the crite- 
rion for state mixing in Cg symmetry. When the symmetry of the Hamiltonian 
is reduced to Sq (C3 upon application of H/ /z) the basis functions may be 
grouped into three different irreducible representations: |0) , |±3) , |±6) e A; 
1+4) , 1+1) , 1-2) , 1-5) e 1^2; 1-4) , 1-1) , 1+2) , 1+5) e Tg. Here, the selection 
rule for mixing is Am^ = 3n, with n being integer, again in agreement with 
the preceding calculations. 

Before concluding this section, the influence of the exchange coupling, J, on 
the QTM observed in Mna deserves further consideration. The J dependence 
of higher-order {p > 4) coefficients in the GSA has been discussed previously 
for several other high-symmetry SMMs [16,18,20,27,28]. In these cases, the 
2'^'^-order transverse anisotropy (q > 0) cancels exactly, emerging at higher 
orders as a consequence of the mixing of spin states. This is illustrated for 
Mn3 in Fig. 1.6, which plots the power law dependence of several QTAf gaps as 
a function of the ratio of J/d; the single-ion zfs parameters given above were 
employed in these calculations. It is found that the QTM gaps are proportional 
to |J|-", i.e., Bf cx |J|-i and cx |J|-2 [18]. Note that this implies a 
complete suppression of QTM in the strong coupling limit (| J| \d\). 

1.2.3 The influence of disorder on QTM 

An important conclusion of the preceding analysis is the demonstration of the 
existence of fc-odd QTM resonances, i.e., a quite realistic parameterization of 
Eqn. 1.6 generates zfs terms in the GSA containing odd powers of >S'+ and S-. 
These ideas should apply quite generally. For example, the disorder potential 
associated with the distortion of a symmetric molecule likely contains zfs 
terms {e.g. 0\ or 0\) that unfreeze fc-odd QTM resonances (as explicitly 
demonstrated in Sec. 1.3.3), contrary to the belief that odd QTM resonances 
cannot be generated in this way [29]. However, it remains to be seen whether 
this can account for the absence of selection rules in SMMs such as Mni2. 
We note that these arguments do not apply to zero-field (fc = 0) QTM in 
half-integer spin systems, which is strictly forbidden according to Kramers' 
theorem [30]. 

This revives a partly unresolved and somewhat controversial issue con- 
cerning the influence of disorder on the QTM characteristics of SMMs. Dis- 
order became a focus of attention in some of the early spectroscopic inves- 
tigations of the Mni2-acetate and FegBr SMMs, revealing significant distri- 
butions (or strains) in the measured GSA D parameters [31 34]. Aroimd the 
same time, Chudnovsky and Garanin argued that long-range strains nucle- 
ated by line dislocations could give rise to a broad distribution of transverse 
2"'^-order anisotropics in otherwise high-symmetry crystals of SMMs such as 
Mni2-acetate, i.e., a broad distribution (on a logarithmic scale) in E centered 
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about an average value of zero [35,36]. Importantly, the Chudnovsky/Garanin 
tlic'ory pointed out that disorder would lead to local variations in molecular 
symmetry away from the ideal {S4 for Mni2-acetate) , and that this could mod- 
ify the selection rules governing QTM. This motivated intense efforts aimed 
at carefully studying QTM in Mni2-acctatc, including selective hole-burning 
experiments targeted at subsets of molecules belonging to different parts of 
the relaxation time distribution [37 40]. A breakthrough was achieved as a 
result of crystallographic investigations by Cornia et al. [41], that revealed 
a form of intrinsic disorder associated with the acetic acid solvent that co- 
crystallizes with the standard form of Mni2-acetate (we refer the reader to 
Refs. [13,39 43] for detailed discussion). The acetic acid forms a hydrogen- 
bond to the Mni2 core, resulting in a non-trivial distortion of the molecule. 
However, while each solvent molecule occupies a position between two Mni2's, 
it can only hydrogen-bond to one of them, with 50:50 probability. Hence, real 
Mni2-acetate crystals contain a statistical distribution of several different sol- 
vent isomers, some of which maintain approximate four-fold symmetry, while 
more than 50% have a lower (rhombic) symmetry [41]. EPR, inelastic neu- 
tron scattering and magnetic hysteresis measurements subsequently yielded 
excellent qualitative and quantitative agreement with the model proposed by 
Cornia, thus demonstrating for the first time that solvent disorder can have 
a profound influence on QTM relaxation [13,39,40,42,43]. 

Many more recent studies have reinforced the idea that solvent disorder 
can significantly infiuence QTM relaxation in SMMs. First of all, magneti- 
zation and EPR studies have shown that the anomalous distributions in zfs 
parameters found for Mni2-acetate are absent in several newer high-symmetry 
(5*4) Mni2 SMMs that do not suffer from the intrinsic solvent disorder (or for 
which the interaction between the solvent and the SMM core is far weaker 
than in the original acetate) [24,28,44-47]. Interestingly, the deliberate re- 
moval of solvent from the newer Mni2's (by pumping on the samples at room 
temperature) has been shown to accelerate the low temperature magnetization 
relaxation, without affecting the height of the classical relaxation barrier [24] . 
Meanwhile, EPR studies demonstrate that the solvent loss induces disorder 
that looks very similar to the intrinsic disorder in Mni2-acetate [24]. This 
again suggests that the induced (extrinsic) disorder causes the faster relax- 
ation, presumably as a result of quantum tunneling processes. This leads to 
known sample handling problems, i.e., crystals containing volatile solvent {e.g. 
Mni2BrAc-CH2Cl2) can change beyond recognition as far as their QTM and 
EPR characteristics are concerned if they arc cooled under vacuum [24,47.48]. 

The reason why the ideal Mni2 SMM is so susceptible to disorder is be- 
cause it has such a high symmetry; the nominally forbidden 2"'*-order trans- 
verse anisotropy rapidly rcemerges upon the introduction of weak disorder, 
either through solvent loss or otherwise. This is not the case for lower sym- 
metry molecules that already possess a 2'"*-order rhombic zfs interaction [49] . 
This has caused some confusion in the literature. As an aside, we note that 
internal transverse dipolar /hyperfine fields can, in principle, also affect QTM 
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selection rules in high-spin SMMs [50] . Indeed, early work demonstrated that 
a combination of allowed transverse zfs interactions, together with transverse 
dipolar/hyperfine fields, may explain the observed absence of QTM selection 
rules in Mni2-acetate in the thermally activated regime [51,52]. However, this 
explanation fails in the pure QTM regime, where tunneling couples low lying 
spin states at a high order of perturbation theory (Amg >> 1). More recent 
studies claim that differences in QTM relaxation observed for the various 
high-symmetry Mni2 SMMs are due entirely to differences in the widths of 
the dipolar field distributions, which obviously depend on the crystal struc- 
tures [50]. However, the recent preparation of two versions of Mni2-acetate 
that are identical in almost all respects (including the lattice constants), apart 
from the co-crystallizing solvent (acetic acid in one case, methanol in the 
other), seem to rule out this assertion [24]. While dipolar fields undoubtedly 
play a crucial role in the collective QTM relaxation in SMM crystals [53] , the 
marked differences in relaxation rates found for the two Mni2-acetates appear 
to be related to the disorder associated with the hydrogen bonding acetic acid 
solvent, which is not present in the methanol variant. 

The large dimension of the Mni2 MS Hamiltonian presents a considerable 
challenge in terms of gaining theoretical insights into the effects of disorder. 
However, there exist many smaller molecules with equally high symmetry 
which are, thus, more amenable to this type of study. Indeed, this issue is 
revisited in Sec. 1.3.3, which deals explicitly with the distortion of a Ni4 SMM 
that possesses the same intrinsic 5*4 symmetry as the ideal Mni2's [16]. Aside 
from the obvious computational advantages, several smaller SMMs are also 
known to crystallize with no lattice solvent molecules [18,22,54]. More impor- 
tantly, there exist families of low-nuclearity SMAIs for which some members 
co-crystallize with solvent, while others do not. These include Mua, Mn4 and 
Ni4, which represent the focus of the remainder of this chapter. The spectro- 
scopic differences between solvatcd and solvent-free SMMs are quite dramatic. 
For instance, D-strain is almost absent in the latter, giving rise to remark- 
ably sharp EPR spectra. This again implicates solvent molecules as a major 
source of disorder in SMM crystals. The key finding involved a solvent-free 
Mn3 compound, which is the only SMM to display a complete absence of 
a symmetry-forbidden QTM resonance [21]. When combined with the ob- 
servation of uniquely sharp EPR spectra [22], this result suggests that it is 
the absence of solvent that unmasks the intrinsic QTM selection rules, again 
hinting at the connection between solvent, disorder and the absence of QTM 
selection rules in other SMMs. 

Fig. 1.7a shows derivatives of magnetization hysteresis curves for Mna, 
recorded at different temperatures from 0.3 K to 2.6 K, with H jj z. At low 
temperatures, the fc = 1 resonance is completely absent. It eventually appears 
for temperatures above 1.3 K as a result of a symmetry allowed thermally 
activated QTM process. As discussed above, the trigonal symmetry of this 
molecule enforces the [Am] = 3n selection rule when taking into accoimt the 
8.5 degree misalignment of the JT axes from the molecular C3 axis {Sq reduced 
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to C3 when a longitudinal field is applied). The effect can be seen in Fig. 1.7b, 
which shows the tunnel splittings for the four lowest resonances /c = — 3, 
calculated by diagonalization of the MS Hamiltonian of Eqn. 1.6 with the pa- 
rameters given in Ref. [21]. In the absence of a transverse field {Ht = 0), the 
ground state tunnel splitting is always absent for resonances k = 1 and k ~ 2, 
while the degeneracy is only broken in resonances A; = and A; = 3. Conse- 
quently, one expects steps in the hysteresis curves (peaks in the derivatives) 
appearing only at fc = and k = Z. The absence of the k = 1 resonance at 
low temperatures constitutes direct evidence for the expected QTM selection 
rule, an observation made possible because of the highly ordered solvent-free 
crystal structure. Following the same reasoning, resonance k = 2 should also 
be absent at low temperatures, since the ground tunnel splitting couples spin 
states differing by |Am| 7^ 3n. However, as can be seen in Fig. 1.7b, the k = 2 
splitting grows quickly as a function of the transverse field, reaching observ- 
able magnitudes for field values provided by internal dipolar fields. This result 
shows that internal Zeeman fields can indeed unfreeze some QTM resonances, 
but not all of them. Indeed, the ground state tunnel splitting associated with 
the k = 1 resonance remains almost two orders of magnitude smaller than 
that of fc = 2 for the same transverse field. One would expect the influence 
of dipolar fields to diminish further still in the pure QTM regime for SMMs 
with larger S values. 

The Mna SMM illustrates perfectly how crystalline disorder can mask the 
fundamental QTM behavior in SMMs; in this particular case, it is the absence 
of disorder that unmasks intrinsic symmetry-enforced quantum properties. 
This, in turn, allows fundamental insights into the influence of the internal 
molecular degrees of freedom on the QTM phenomenon. The low-nuclcarity of 
the Mn3 SMM proved particularly helpful by making this a computationally 
tractable problem. The following section digs deeper into the unusual BPI 
patterns predicted for trigonal SMMs. 

1.2.4 Berry Phase Interference in Trigonal Symmetry 

This section focuses explicitly on the BPI patterns generated by the Of op- 
erator. In contrast to all of the evcn-q GSA terms, the xy-plane does not 
correspond to a sj^mmctry clement for the odd-(jr interactions, as discussed in 
Sec. 1.1.1. Hence, the O4 operator is expected to result in BPI patterns which 
have not been observed in previous studies of SMMs, essentially all of which 
possess even rotational symmetry with only even-g zfs interactions [55-57]. 

The influence of S4O4 is quite fascinating. In order to simplify discussion. 
Fig. 1.8 was generated with Bg ~ 0; details of the calculations can be found 
in Ref [58]. The Ag g (fc = 0) QTM resonance exhibits the most intriguing 
new features. One might expect a six-fold behavior due to the intrinsic Cj 
symmetry of the Hamiltonian, i.e., the spectrum should be invariant under 
inversion of iJx- However, this assumes that Hj^ = 0. In fact, application of 
a transverse field causes a shift of the Ag^g minimum away from Hj, = 0, as 
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illustrated in Fig. 1.8b, i.e., the k = QTM resonance shifts away from Hi, = 
in the presence of a finite transverse field. The resultant transverse-field BPI 
patterns appear to exhibit a hexagonal symmetry in Fig. 1.8a. However, the 
color coding represents the polarity of the required compensating longitudinal 
field, Hl- Thus, on the basis of the sign of Hi,, one sees that the BPI minima 
in fact exhibit a three-fold rotational symmetry, which is consistent with the 
symmetry of the BlOl interaction. One way to interpret this result is to 
view the 0| operator as generating an effective internal longitudinal field, 
-ff£, under the action of an applied transverse field; is then responsible 
for the shift of the k = resonance from i?L = 0. This can be seen from 
the expression of the Of = ^[Sz,S^ + S^] operator, which, unlike the g-cvcn 
operators, contains an odd power of Sz, akin to the Zeeman interaction with 
i? II 2. An alternative view may be derived from the Sq surface depicted in 
Fig. 1.2c, where one sees that the hard/medium directions do not lie within 
the a;y-plane, contrary to the case for the g-even operators. In other words, 
the classical hard plane is not flat, but corrugated with a 120° periodicity. 
Consequently, application of a longitudinal field is required in order to insure 
that the total applied field lies within the hard plane when rotating H^. 
Note that the predicted BPI patterns nevertheless exhibit the required Cj 
symmetry, i.e., they arc invariant with respect to inversion of the total field. 

Fig. 1.8b plots the shift of the A; = resonance (Ag g minimum) away 
from Hi, = upon applying a transverse field, Ht, for several orientations 
within the xy-plane. The shift is positive for 0° and 15°, and negative for 
45° and 60°, with no shift at 30° {i.e. the 30° resonance occurs at Hj, = 0). 
In other words, the quantum molecular hard plane is not flat, but rather 
corrugated, with a 120° periodicity. This is consistent with the classical energy 
surface shown in Fig. 1.2c. It is also notable that the Hi, shift displays a non- 
linear dependence on Ht, which indicates that the exact locations of the hard 
directions depend on the magnitude of H^. Finally, it should be emphasized 
that these phenomena, especially the shift of the fc = resonance from Hi, = 0, 
cannot be generated by any of the even-g operators [57] , where the xy plane 
necessarily corresponds to the hard plane because of the additional mirror 
symmetry about this plane (see the discussion in Sec. 1.1.1). 

1.3 Quantum Tunneling of Magnetization in the 
High-Symmetry Ni4 Single-Molecule Magnet 

1.3.1 The Ni4 Single-Molecule Magnet 

The [Ni(hmp)(dmb)Cl]4 SMM (henceforth Ni4) possesses q-even rotational 
symmetry [14, 16,54,59,60]. The complex crystallizes in an /41/a space group 
without any lattice solvent molecules. The structure of the Ni4 molecule is 
shown in Fig. 1.9a. The magnetic Ni404 core is a slightly distorted cube with 
the Ni^^ ions (s = 1) located on opposite corners, as sketched in Fig. 1.9b. The 
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distorted cube retains ^4 symmetry, with the 54-axis indicated in Fig. 1.9a. 
The four Ni^^ ions arc fcrromagnctically coupled, leading to a spin S* = 4 
molecular ground state. The Ni4 SMM exhibits extremely fast zero-field QTM, 
which significantly reduces the effective relaxation barrier. Nevertheless, it 
does display a small magnetic hysteresis [60]. However, the fast relaxation un- 
fortunately precludes the observation of fc > QTM resonances. Nevertheless, 
a theoretical study of Ni4 proves enlightening. The molecule can be described 
with the following spin Hamiltonian: 



This Hamiltonian differs from Eqn. 1.6 in that the individual rotation ma- 



trices, Ri, are replaced by a single matrix, n, which explicitly takes the ro- 
tational symmetry of the molecule into account, including cases involving 
improper rotations. The zero-field anisotropy can then be parameterized by 

a single rf-tensor (corresponding to one of the ions), specified with respect 

to the molecular coordinate frame. If the local coordinate frames of the in- 
dividual ions are related by a series of proper rotations (C2, C4, etc.) within 
the molecular coordinate frame, then n may be replaced by a single rotation 

matrix R (corresponding, e.g., to a 90° rotation about z for a molecule with 
C4 symmetry). On the other hand, if the local coordinate frame of the i^^ ion 
is related to the molecular frame by an improper rotation (5*4 or Sf), then 

n= a R, where a represents a reflection in the plane perpendicular to the ^4 

2 4 
axis. Note that, for S4 symmetry, n is equivalent to a C2 rotation, and n is 

equivalent to the identity matrix. 

The Ni4 SMM is a particularly ideal platform for comparison with Mns. 

The molecule possesses a well separated 5 = 4 ground state with the 5 = 3 
excited spin multiplets located roughly 30 K above in energy. The 3x3 Hamil- 
tonian matrix associated with a single Ni^' ion contains only two 2"'^-order 
zfs parameters, d and e, i.e., higher order single-ion anisotropics (p > 4) arc 
strictly forbidden. The zfs of the individual Ni^^ ions, as well as their orienta- 
tions, have been directly measured through EPR studies on an isostructural 
diluted Zn4_,-,.Ni;,. compound (sec Fig. 1.10(a) and Ref. [19]). Due to the restric- 
tion of 5*4 symmetry, only two independent Heisenberg interaction parameters, 
Ji and J2, are allowed; these interactions can be determined by dc susceptibil- 
ity measurements [60]. Therefore, all of the parameters in Eqn. 1.9 are known 
independently. Meanwhile, the molecule possesses the same 54 symmetry as 
Mni2, which prohibits the rhombic anisotropy term in the GSA Hamilto- 
nian. The high symmetry of the molecule has been confirmed by single-crystal 
EPR measurements, where exceptionally sharp resonances are again observed, 
with a four-fold modulation pattern upon rotating Ht (see Fig. 1.10(b) and 
Ref. [61]). This clearly illustrates the presence of a 4*^- (or higher-) order 
transverse GSA anisotropy which is responsible for the fast QTM. 



4 
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1.3.2 Quantum Tunneling of Magnetization in the Ni4 SMM 

In analogy with Mn.3, the transverse GSA anisotropy in Ni4 is assessed by cal- 
culating the QTM gaps, focusing on the fc = 0, 1 . . . 4 ground state resonances, 
as shown in Fig. 1.11; the |Am| values associated with these resonances equal 
8, 7 ... 4, respectively. The simulations were performed using the published zfs 
parameters d = -7.6 K, e = 1.73 K and Ji = J2 = -10 K [14,16,19]. Previous 
EPR studies also show that the local easy-axes of the Ni^^ ions arc tilted away 
from the molecular ^-axis by 6 = 15° (See Fig. 1.9c). However, the 9 = 0° case 
is also examined in order to further explore the influence of easy-axis tilting 
on the symmetry of the molecular Hamiltonian. 

Fig. 1.12 shows the ground state QTM gaps as function of transverse field 
(Ht), deduced via exact diagonalization of Eqn. 1.9. As seen in the figure, 
^4.4 {f^- = 0) and A0.4 (fc = 4) retain non-zero values in the absence of a 
transverse field, while all other tunnel splittings vanish at Ht = 0. This result 
is not surprising based on the 54 molecular symmetry, where only |Am| = in 
(n is an integer) QTM resonances are allowed. However, unlike the Mn^ SMM, 
the QTM selection rules corresponding to the 9 = 15° and 0° situations are 
exactly the same. In both scenarios, only the Aj 4 {k = 0) and Ao,4 {k = 4) 
gaps are non-zero, while the other fc-even QTM gap, A2,4 {k = 2), vanishes 
when Ht = 0. These results imply that the easy-axis tilting does not affect the 
symmetry of the Hamiltonian, contrary to the case for the Mn^ SMM. This 
can be understood in terms of the different symmetry properties associated 
with g-even and q-odd cases. In the even case, the molecular z-axis must also 
be a C2 axis. Consequently, forcing the local C2 axes of the individual ions 
to be parallel to the molecular z-axis {9 = 0°) does not introduce an extra 
C'2 symmetry to the molecular Hamiltonian. In contrast, the molecular z-axis 
is not a C2 axis in a molecule with odd rotational symmetry. Therefore, the 
symmetry of the molecular Hamiltonian changes when 9 = 0°. 

In the preceding discussions of Mns , the QTM selection rules can be sim- 
ply understood in terms of the rotational symmetry of the molecule (Cg or 
C3). In contrast, the selection rules for Ni4 cannot be fully explained by the 
5*4 molecular symmetry; one must additionally take into account the intrinsic 
Ci symmetry of the spin Hamiltonian. Upon application of a magnetic field 
parallel to the molecular z-axis, the 6*4 symmetry group reduces to C'2, for 
which the A2_4 (k = 2) QTM resonance should be allowed. This clearly contra- 
dicts the simulation in Fig. 1.12, which suggests a higher symmetry. However, 
one must also consider the Ci symmetry associated with the SO interaction. 
The consequential zero-field spin Hamiltonian then possesses >S'4 x Cj = C^h 
symmetry, which corresponds to the symmetry of the Of interaction, as seen 
in Fig. 1.2d. Upon application of a longitudinal field, the C4/1 group reduces 
to the C4 group, for which the expected QTM selection rule |Am| = 4n is 
recovered. The Ci symmetry is guaranteed by the nature of the SO inter- 
action. This property is not limited to spin Hamiltonians, i.e., it applies to 
any Hamiltonian dictated by crystal field and/or SO physics, where the Q 
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symmetry should apply regardless of whether the orbital angular momentum 
is qucnichcd or not. In other words, it is always necessary to consider the Ci 
symmetry in addition to the structural symmetry, especially when improper 
rotations are involved. Unfortunately, observation of A; > QTM steps in Ni4 
is impractical due to the extremely fast tunneling at fc = 0. This tunneling 
should be greatly suppressed if the ground spin state of the molecule is in- 
creased, as is of course the case for Mni2. However, it would be interesting 
to obtain a four-fold symmetric SMM constituted of four s = 2 Mn^^^ ions, 
for which it would be possible to study the A; > QTM steps. Moreover, the 
Hamiltonian dimension of just 625 x 625 would be quite manageable. 

1.3.3 Disorder 

In the presence of random disorder, one would expect the symmetry of most 
molecules to be lowered, leading to an absence of QTM selection rules. The 
Ni4 molecule provides an excellent platform to study this issue. Fig. 1.13 was 
generated by adjusting the orientation of the zfs tensor of one of the Ni^^ ions in 
the molecule, i.e., the zfs tensors of three of the Ni ions are tilted 15° from the 
molecular z-axis, while the other is tilted 10°. It should be emphasized that it 
is not trivial to find the orientation of the molecular easy-axis in this situation, 
i.e., it no longer coincides with the molecular 2;-axis. For each resonance {k = 
to 4), a search was performed for the minimum QTM gap by varying the 
orientation of the applied field. As seen in the figure, all resonances posses a 
non-zero QTM gap at Ht = 0. The inset to Fig. 1.13 plots A3_4 (A; = 1), A2,4 
(fc = 2) and Aj 4 {k = 3) on a log- log scale, clearly demonstrating that these 
QTM gaps, which are forbidden for ^4 symmetry, also saturate at non-zero 
values when Ht — ^ 0. These results show that a small disorder effectively 
unfreezes all QTM steps without the assistance of a transverse field. This 
argument can be reinforced by group theoretic considerations. With random 
disorder, the symmetry of a molecule is lowered to Ci, resulting in a spin 
Hamiltonian with Q symmetry. Upon applying a longitudinal field, the Ci 
group reduces to Ci, where all of the states necessarily belong to the same 
one-dimensional irreducible representation [25]. Therefore, mixing between all 
states is allowed. We note that this kind of disorder can be introduced by small 
crystallographic defects, which always exist to some degree in real samples. 
Thus, exceptionally clean crystals are required in order to observe symmetry 
imposed QTM selection rules. Importantly, the preceding discussion clearly 
demonstrates that disorder can be responsible for the observation of A;-odd 
QTM steps in SMMs with even rotational symmetries. 
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1.4 Quantum Tunneling of Magnetization in 
Low-Symmetry Mn4 Single-Molecule Magnets 

In order to contrast results presented in previous sections, EPR and QTM/BPI 
results are presented here for two related Mn4 SMMs that possess almost no 
symmetry. Both molecules crystallize in the triclinic PI space group. One of 
the structures co-crystallizes with solvent, while the other does not. Conse- 
quently, significant differences are observed in terms of the widths of EPR and 
QTM resonances due to the different degrees of disorder in the two crystals. 
In addition, small structural differences associated with the Mn4 cores result 
in different coupling strengths between the Mn ions within the two molecules 
which, in turn, result in different QTM behavior. 

1.4.1 The Mn4 Single-Molecule Magnets 

The Mn4 molecules (Figs. 1.14a and b) possess mixed- valent butterfly- type 
structures, with two central Mn^^' ions (s2 = S4 = 2) in the body positions and 
two Mn'^ ions (,si = S3 = 5/2) in the wing positions (see Fig. 1.14c). Magnetic 
superexchange is mediated through oxygen bridges. The Mn4 molecule that 
co-crystallizes with solvent is [Mn4(anca)4(Hedea)2(edea)2]-2CHCl3-2EtOII, 
henceforth Mn4-anca (see Fig. 1.14a and Ref. [8] for more details). The solvent- 
free molecule is [Mn4(Bet)4(mdea)2(mdeaH)2](BPh4)4, henceforth Mn4-Bet 
(see Fig. 1.14b and Ref. [9] for more details). Both molecules crystallize in 
the ccntrosymmetric triclinic space group PI. The asymmetric unit therefore 
consists of half the molecule (Mn^^^Mn^^), with the other half generated via an 
inversion. This also ensures that all four Mn ions lie in a plane. Temperature 
dependent susceptibility data recorded at different magnetic fields indicate 
that both molecules possess a spin S = 9 ground state at low temperatures, 
implying overall ferromagnetic coupling within the molecules (note that this 
does not rule out the possibility that one of the exchange paths could be 
antiferromagnetic) . 

1.4.2 EPR and QTM Spectroscopy in Mn4 SMMs with and 
without Solvent 

A selection of EPR and QTM measurements from Refs. [8, 9] are presented 
in Fig. 1.15: (a) displays 165 GHz EPR spectra recorded at different tempera- 
tures for a single crystal of Mn4-anca, with the magnetic field applied close to 
the molecular easy-axis; (b) shows equivalent spectra obtained for Mn4-Bct 
at a frequency of 139.5 GHz and similar temperatures. The first thing to note 
are the obvious differences in the EPR linewidths in the two figures. This 
again provides a clear illustration of the inferior quality of samples in which 
the SMMs co-crystallize with interstitial solvent molecules. In the present ex- 
ample, the Mn4-anca sample is the more disordered, resulting in a broader 
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distribution of GSA zfs parameters. A series of nine absorption peaks can 
clearly be seen for the Mn4-anca SMM in Fig. 1.15a, which have been labeled 
1 through 9. These correspond to transitions from consecutive spin projection 
(ms) states belonging to the 5 = 9 ground state multiplet, where the num- 
bering denotes the absolute ms value associated with the level from which 
the EPR transition was excited, e.g., resonance a = 9 corresponds to the 
ms = — 9 — >■ — 8 transition. The fact that all of the spectral weight transfers 
to the a = 9 resonance as T — > indicates uniaxial anisotropy, i.e., D < 
according to the GSA Hamiltonian of Eqn. 1.2. The uneven spacing between 
the labeled ground state resonances is indicative of 4**^- (and higher-) order 
anisotropy within the GSA (or weak exchange within the MS picture). Finally, 
as the temperature is increased, a few weaker resonances (not labeled) can be 
seen to appear in Fig. 1.15a between the labeled transitions. These additional 
peaks are associated with the population of higher-lying, excited spin states, 
e.g., S = S. 

The EPR spectra obtained for Mn4-Bet (Fig. 1.15b) are not so simple to 

interpret. First and foremost, many more peaks are observed, suggesting the 
population of many more spin states. Based upon the knowledge gained from 
Mn4-anca, and the results of subsequent simulations, the nine resonances cor- 
responding to transitions within the S = 9 groimd state multiplet have been 
labeled in the figure. Meanwhile, the peaks that are not labeled correspond 
to transitions within low- lying excited spin multiplets. Clearly, the emergence 
of excited state EPR transitions at much lower temperatures indicates signif- 
icantly weaker exchange coupling in the Mn4-Bet molecule; note that excited 
state intensity is seen even at the lowest temperature, whereas this is not the 
case until ^ 9 K in the Mn4-anca sample. The weaker exchange and higher 
crystal quality associated with the Mn4-Bet sample lead to the observation of 
unusual MQT/BPI behavior, as detailed in the following section. 

Although one can reproduce the peak positions of the nine labeled EPR 
transitions in both Figs. 1.15a and b using the GSA (including p > 4 terms), 
a MS Hamiltonian becomes essential to account for transitions within ex- 
cited spin states. In other words, one starts to see the limitations of the GSA 
in these two cases— especially for Mn4-Bet. Diagonalization of the exact MS 
Hamiltonian that considers all four Mn ions and the couplings between them 
(as indicated in Fig. 1.14c) is manageable on a standard computer. However, a 
convenient and reasonable approximation replaces the ferromagnetically cou- 
pled central Mn''^ ions with a single Ss = 4 spin, resulting in a linear trimer 
consisting of the central spin, sb, and the two outer = 5/2 Mn^^ spins, as 
depicted in Fig. 1.14d. This approximation, which contains elements of both 
the GSA and MS models, has been compared with the more exact four-spin 
model for the case of Mn4-Bet in Rcf . [7] . The two models give good agreement 
in terms of the lowest-lying portions of the energy-level diagrams that dictate 
the low-temperature QTM and EPR properties, provided that the ferromag- 
netic coupling between the Mn^^^ ions is not too weak. Although it involves a 
level of approximation, the trimer model captures all of the important physics 
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associated with these low-symmetry SMMs. Moreover, the smaller Hamilto- 
nian matrix dimension enables much faster fitting (hours instead of days) , and 
employs fewer parameters— just two d tensors and a single exchange coupling 
constant, J. Finally, working with a single J parameter to identify the internal 
exchange energy becomes particularly useful in relating the distinct behavior 
of the two molecules to different intramolecular couplings which result from 
slight structural disparities between the two compounds. In the following, we 
diagonalize the MS Hamiltonian (Eqn. 1.6) using the trimcr model (Fig. 1.14c) 
to account for the energy landscape associated to the lowest lying molecular 
spin multiplets, which result from the main anisotropy terms in the Hamilto- 
nian (i.e. axial terms). The full MS Hamiltonian (Eqn. 1.6) including the four 
manganese ions (Fig. 1.14d) is used to account for the behavior of the tunnel 
splittings, which result from the smaller anisotropy terms in the Hamiltonian 
(i.e. transverse terms) and are more sensitive to small variations of the internal 
degrees of freedom of the molecules. 

In fitting the data in Figs. 1.15a and 1.15b, as well as other EPR data ob- 
tained at different temperatures and applied field orientations (see Refs. [8,9]), 
one finds that the exchange coupling constant J has a strong influence on the 
positions of the EPR peaks (particularly the relative spacings between peaks). 
This again highlights the fact that one cannot use a GSA to realistically de- 
scribe these results, especially those of Mn4-Bet, i.e., there is no exchange 
parameter in the GSA (all energy splittings are parameterized in terms of the 
SO anisotropy). It is thus preferable to use the MS approach whenever compu- 
tational resources will allow, as is the case for all of the low-nuclearity SMMs 
described in this chapter. Indeed, there are many other interesting features 
associated with the magnetic behavior of SMMs that cannot be explained 
with the GSA, regardless of how complex the corresponding Hamiltonian is, 
e.g., QTM resonances involving level crossings between different spin states. 

QTM spectroscopy also facilitates comparisons between the two Mn4 
molecules. Figs. 1.15c and d show the association between the observed QTM 
resonance positions and the corresponding energy-level diagrams for the Mn4- 
anca and Mn4-Bet molecules, respectively. The QTM resonances are deter- 
mined from the positions of the peaks in the derivatives of the magnetization 
versus field curves obtained at different temperatures, as shown in the lower 
portions of Figs. 1.15c and d. Note that the effects of the solvent disorder can 
again be seen, causing broader QTM resonances for Mn4-anca in comparison 
to Mn4-Bet. The energy level diagrams are obtained via exact diagonaliza- 
tion of the MS Hamiltonian of Eqn. 1.6 using the trimer model depicted in 
Fig. 1.14c with the following parameters: J = 5.42 K, di = d2 = cLa = —0.115 
K, ds = -2.22 K in Mn4-anca (Fig. 1.15c); and J = 1.90 K, = -0.115 
K, ds = 2.00 K in Mn4-Bet (Fig. 1.15d), with isotropic g = 2.0 for all ions 
in both cases. The correspondence between QTM resonances and their as- 
sociated level crossings are indicated by vertical arrows. Crossings involving 
the ground state jS* = 9, = 9) with other levels within the same multiplet 
{S = 9) are indicated by black arrows and expected to appear at the lowest 
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temperatures, for which transitions involving all excited states (red arrows) 
should vanish. Interestingly, some of the resonances (blue arrows) correspond 
to crossings between levels associated with different spin states, i.e., S = 9 
and 5 = 8. The main difference between the two molecules resides in the value 
of J, being more than double for Mn4-anca. This results in the lowest spin 
projection states (mg = ±8) associated with the first excited state {S = 8) 
being much closer to the ms = ±9 ground states in Mn4-Bet (~ 8 K sepa- 
ration) than in Mn4-anca (^ 22 K separation). These findings are consistent 
with the temperature dependence of the EPR spectra, which suggested the 
excited states to be appreciably lower in energy in the Mn4-Bet molecule in 
comparison to Mn4-anca. The differences in J values can be reconciled with 
the minor structural differences between the two molecules. It is well known 
that the superexchange coupling between two transition metal ions is very 
sensitive to the bridging angle, to the extent that the sign of the interaction 
can switch from ferromagnetic to antiferromagnetic within a small range of 
angles [18,62,63]. Indeed, there are measurable differences in the bond angles 
and distances associated with these two Mn4 molecules. 

1.4.3 Berry Phase Interference in Mn4-Bet 

The spectroscopic results presented in the previous section illustrate how small 
structural perturbations can lead to significant changes in the exchange cou- 
pling within a molecule. Crucially, the Mn4-Bet SMM resides in a particularly 
interesting region of the 'anisotropy' versus 'exchange' parameter space in 
which excited spin states exert a significant influence on the low-energy /low- 
temperature quantum properties. First and foremost, it can be seen that some 
of the QTM resonances involve level crossings between different spin multi- 
plets. More importantly, the QTM properties within this intermediate ex- 
change regime {d ~ J) are extraordinarily sensitive to the internal magnetic 
structure of the molecule. As noted in previous sections, the physics associated 
with the strong exchange limit (J >> d) is dominated by the 2"^-order GSA 
anisotropy. Consequently, any observable BPI patterns should display a high 
degree of symmetry {D2h), regardless of the molecular symmetry. However, in 
the intermediate exchange regime, one may expect any BPI effects to mimic 
the symmetry of the molecule under investigation much more closely. This 
symmetry can be expressed exactly using a four spin MS model in the case of 
Mn4-Bet, although we note that one can also reproduce some of the features 
of these experiments using a GSA that includes appropriate 4"^- (and higher-) 
order terms. The virtual lack of any symmetry associated with Mn4-Bet leads 
to some remarkable BPI patterns, which also shed light on a previous mystery 
surrounding another high-nuclearity Mni2 SMM [64]. 

Low temperature QTM measurements performed in the presence of a 
transverse magnetic field, H^, enable exploration of the dominant symme- 
tries associated with the Hamiltonian describing a SMM. Fig. 1.16a shows the 
modulation of the QTM probability for the fc = resonance for Mn4-Bet as 
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a function of the magnitude of Ht applied parallel to the hard axis of the 
molecule ((/) = 0°, sec Ref. [7] for details). As explained in Sec. 1.2.4, the os- 
cillations correspond to BPI (constructive/destructive interference associated 
with equivalent tunneling trajectories on the Bloch sphere), with minima oc- 
curring at regularly spaced field values (Affx = 0.3 T). Experiments designed 
to probe the modulation of the fc = QTM gap as a function of the orientation 
of a fixed transverse field within the hard plane (see Ref. [7]) reveal a two- 
fold behavior. One may be tempted to ascribe this to a rhombic anisotropy. 
However, the molecule possesses a much lower symmetry (PI). In fact, the 
two-fold pattern is a direct manifestation of the Cj symmetry associated with 
the SO interaction. Since no longitudinal field [Hl) is present for the fc = 
resonance, the Hamiltonian must be invariant with respect to inversion of 
i?T— hence the apparent two-fold behavior. Note that the fc = BPI oscil- 
lations do, indeed, respect the symmetry of the Hamiltonian, i.e., they are 
invariant under inversion of Ht. 

Due to the absence of Hi,, fc = turns out not to be the most inter- 
esting QTM resonance, because the Ci symmetry guarantees symmetric BPI 
patterns about ifx = 0. In contrast, this is clearly not the case for the BPI 
pattern associated with the fc = 1 resonance, as can be seen in Fig. 1.16b. In 
this case, a single interference minimum is observed at iJx = 0.3 T for only 
one polarity of the transverse field, i.e., the corresponding BPI minimum is 
completely absent under inversion of H^. Such a result is not so surprising 
when one recognizes that there are no mirror symmetries within PI. Hence, 
there is no reason why the BPI patterns should be invariant under inversion 
of just one component of the applied field. However, the Hamiltonian, and 
therefore the BPI patterns, must be invariant under a full inversion of the 
applied field, i.e., inversion of both ffx and Hi, together. This indeed turns 
out to be the case for the fc = 1 resonance, as clearly seen in Fig. 1.16b. 

Another interesting feature observed in the BPI patterns of Mn4-Bet is 
the motion of the minima associated with the fc = 1 resonance within the 
xy-plane. This can be observed in Fig. 1.16d, which shows a color contour 
plot of the QTM probability in fc = 1, as a function of the magnitude and the 
orientation of ffx- Two minima can clearly be observed; they are again spaced 
by ~ 0.3 T, and are located half way between the fc = minima, as is usually 
the case for even/odd resonances. However, the fc = 1 minima do not appear 
at the same orientations within the xy-plane as those of fc = 0. Moreover, the 
orientations of the two observed fc = 1 minima do not even coincide: (j) = 13.5° 
for the first minimum and <f> = 6° for the second one. Note that, in contrast 
to fc = 1, all of the fc = minima lie along the nominal hard anisotropy 
axis (0 = 0), as seen in Fig. 1.16c. In essence, the hard directions associated 
with the fc = 1 resonance (for which both Hi, and iJx are finite) do not occur 
along a fixed axis, as would be the case for a rhombic molecule. This property, 
which is analogous to the behavior seen in Fig. 1.7b, is a direct consequence of 
the absence of any mirror symmetries in the PI space group. It is impossible 
to simultaneously satisfy both Ci symmetry and a mirror symmetry if the 
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BPI minima do not reside on a fixed axis in space. However, if the mirror 
symmetry is broken, then the BPI minima may in principle occur anywhere 
on the Bloch sphere, so long as they occur in pairs related by inversion. The 
results displayed in Fig. 1.16d were obtained by rotating Ht while keeping 
i?L fixed; the space above and below the xy-plane was not explored due to 
experimental constraints. One cannot rule out further minima above or below 
the xy-plane. Indeed, this may explain why the 2"** k = 1 minimum is so 
weak, i.e., its real location may be above or below the .Ty-planc. In fact, this 
could also be true for the k = minima, as again emphasized in Fig. 1.7b for 
the case of the Mn^ molecule, which lacks xy mirror symmetry. More detailed 
(and time consuming) experiments are clearly required to further explore this 
issue in Mn4-Bet. 

In order to simulate the observed BPI patterns, one must obviously break 

some or all of the mirror symmetries within the Mn4-Bet spin Hamiltonian, 
whilst also respecting the inversion symmetry of the real molecule. There 
really is only one way to achieve this, involving misalignment of the 2"*^- 
order zfs tensors associated with the Mn'^^^ and the Mn^^ ions. The molecular 
inversion symmetry guarantees that the JT axes associated with the Mn^^^ ions 
be parallel to each other; likewise the zfs tensors associated with the Mn^^ ions 
must be co- linear. However, there is no requirement that the tensors associated 
with the two types of ion be co-linear. Indeed, all of the results in Fig. 1.16 
have been reproduced in Ref. [7] following exactly this approach. Although 
the trimcr model (Fig. 1.14c) can reproduce the observed behavior quite well, 
the four spin Hamiltonian (Eqn. 1.6) was employed in order to describe the 
geometry in Fig. 1.14d, since it gives a better quantitative agreement and 
allows for a more physical interpretation of the observations (e.g. the real 
dipolar coupling between the four Mn ions can be employed, which involves 
no fitting parameters). Using this approach, we find the optimal parameter set 
to be as follows: (central Mn^^'s) d2 = = -4.99 K and e2 = = 0.82 K, 
with the easy and hard axes along z (a2 = 0) and x (/32 = 0), respectively; 
(outer Mn'^'s) ^1=^3 = —0.67 K and ei = 63 = 0, with the axes rotated with 
respect to the central spin by identical Euler angles ai_3 = 45°, /3i^3 = 0° (as 
required by inversion symmetry); 7 being zero for all ions; finally, isotropic 
ferromagnetic exchange constants = —3.84 K, = —1.20 K and Jc = 
—3.36 K were employed. It should be stressed that these parameters were 
additionally constrained by the EPR and QTM data in Fig. 1.15. They do not 
necessarily constitute the correct parameterization, but they account for all 
experimental observations. Note that the zfs tensors of the two Mn^^ ions are 
tilted by 45 degrees with respect to the Mn^^^ tensors, thereby breaking the 
molecular xy mirror symmetry. This results in a breaking of the xy mirror 
symmetry of the corresponding spin Hamiltonian. Note that this would not 
be the case for a molecule with even rotational symmetry, because of the Ci 
symmetry associated with the SO interaction. However, in PI {q = 1), the xy 
mirror symmetry of the Hamiltonian is broken, as was also the case for the 
trigonal (g = 3) Mn3 SMM. One could, in principle, also explain these results 
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using a GSA by introducing 4*^- (and higher-) order terms. However, a more 
natural and satisfying account of the results is obtained by diagonalizing the 
four-spin Hamiltonian, which easily allows for a tilting of the zfs tensors of 
the four Mn ions. 

Interestingly, there again exists a connection to Mni2, albeit a wheel 
molecule that bears no resemblance to the well studied, high-symmetry Mni2's 
discussed in Sec. 1.2.3. The Mni2 wheel molecule [64], which possesses the 
exact same PI symmetry as the Mn4 molecules considered in this section, 
attracted considerable controversy on account of the observation of asym- 
metric fc > BPI patterns [17,65]. Initial attempts to account for this be- 
havior involved treating the molecule as a dimer, including an unphysical 
Dzyaloshinskii-Moriya coupling between the two halves of the dimer (this in- 
teraction is forbidden on account of the molecular inversion symmetry [66]). 
The present studies have shown that asymmetric > BPI patterns are, 
in fact, quite natural in low-symmetry molecules. As in the case of Mna, a 
detailed understanding of the QTM characteristics in the Mn4 molecules is 
made possible due to the low nuclearity of the system, which enables the em- 
ployment of a MS Hamiltonian. This, in turn, provides fundamental insights 
that are much harder to achieve when studying larger molecules using a GSA. 

1.5 Summary and Outlook 

This chapter focuses on the microscopic factors that dictate the QTM behav- 
ior observed in polynuclear transition-metal SMMs, with particular focus on 
molecular symmetry. The examples provided involve relatively simple, low- 
nuclearity clusters (Mns, Mn4 and Ni4) which display essentially the same 
physics as the original Mni2 and Peg SMMs that have occupied chemists and 
physicists working in this field for nearly two decades. The simpler systems 
are amenable to analysis using a microscopic spin Hamiltonian that incorpo- 
rates both the single-ion physics, and isotropic exchange coupling between the 
constituent ions, and relies on relatively few parameters. One can therefore 
systematically investigate the role of internal spin degrees of freedom within 
a molecule, in contrast to the approximate GSA approach employed for most 
studies of Mni2 and Fcs. Comparisons between theory and experiment are 
presented for a range of cluster symmetries, with remarkable quantitative 
agreement achieved in all cases. 

Until fairly recently, most SMM research was directed towards polynu- 
clear 3d transition metal clusters, with the synthetic goal of maximizing both 
the molecular spin state and the magneto-anisotropy. However, a number of 
fundamental factors have limited progress based on this strategy, with the 
record blocking temperature for a Mng cluster [67] only just surpassing that 
of the original Mni2 SMM [24]. Limiting factors include: (i) a tendency for 
superexchange interactions between constituent transition metal spins to be 
both weak (few cm~^) and often antiferromagnetic; (ii) the fact that strong 
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crystal-field effects typically quench the orbital momentum associated with 
3d elements, thus significantly suppressing the magncto-anisotropy: and (iii) 
the difficulties associated with maximally projecting any remaining (2"'* order 
SO) anisotropy onto the ground spin state. In fact, careful studies of this issue 
suggest that one is unlikely to achieve anisotropy barriers that significantly 
exceed those of the constituent ions [18]. This is perhaps best illustrated by 
the optimum Mn^, Mhq and Mni2 SMMs, which possess similar barriers (to 
within a factor of < 2 [18]). This is because the molecular anisotropy, D, is 
given by a weighted sum of the anisotropics of the constituent ions (dj), where 
the weighting is inversely proportional to the total molecular spin, S [18,68]. 
Thus, D decreases as S increases, and the theoretical best that one can hope 
to achieve is an anisotropy barrier (~ DS"^) that scales linearly with S (or N, 
the number of magnetic ions in a ferromagnetic molecule). Even then, many 
challenges remain— some fundamental (quantum tunneling, spin state mixing, 
etc. [6]), some synthetic. The synthetic challenges, in particular, become more 
complex with increasing molecule size, e.g., maintaining ferromagnetic cou- 
plings, maintaining parallel arrangements of the individual anisotropy tensors, 
etc. Thus, it is perhaps no surprise that the optimum [Mn^^^jjv SMM has a 
nuclearity of just six [67]! 

Given the aforementioned situation, it has become clear that a more di- 
rect route to magnetic molecules that might one day be used in practical 
devices involves the use of ions that exhibit considerably stronger magnetic 
anisotropics than those that have traditionally been employed in the synthesis 
of large polynuclear clusters, i.e., ions for which the orbital momentum is not 
quenched, and/or heavier elements in which strong SO effects are manifest. 
Examples include certain high-symmetry and low-coordinate 3d transition 
metal complexes (Fe^^ [69], Co^^ [70], even Ni^^ [71]), as well as elements further 
down the periodic table such as the 4/ and 5/ elements. Indeed, over the past 
few years, a number of mononuclear complexes have been shown to exhibit 
magnetization blocking of pure molecular origin [69,70,72-74]. However, the 
quantum magnetization dynamics of these single-ion molecular nanomagnets 
has yet to be studied in detail, and much remains to be learned theoretically 
Obviously, much of the physics related to exchange which is discussed in this 
chapter does not apply in these cases. Nevertheless, the spin Hamiltonian for- 
malism remains applicable, as does the crucial importance of molecular and 
crystallographic symmetry. In particular, 4*'' and higher-order crystal-field in- 
teractions may be expected to play a crucial role in the quantum dynamics 
of mononuclear lanthanidc SMMs [75]. Thus, one would expect similar com- 
binations of theory and spectroscopy to contribute in future to this evolving 
field of research. 



26 Junjie Liu, Enrique del Barco, and Stephen Hill 

I. 6 Acknowledgements 

This work was supported by the US National Science Foundation, grant num- 
bers DMR0804408 (SH), CHE0924374 (SH), and DMR0747587 (EdB). Work 
performed at the National High Magnetic Field Laboratory is supported by 
the National Science Foundation (grant number DMRl 157490), the State of 
Florida and the Department of Energy. 

References 

1. Rudowicz, C. and Misra, S. K. (2001) Applied Spectroscopy Reviews 36(1), 11- 
63. 

2. Gatteschi, D., Sessoli, R., and Villain, J. (2006) Molecular Nanomagnets, Oxford 
University Press, Oxford. 

3. Abragam, A. and Bleaney, B. (1986) Electron paramagnetic resonance of tran- 
sition ions, Dover, New York. 

4. Rudowicz, C. and Chung, C. Y. (2004) Journal of Physics: Condensed Matter 
16(32), 5825. 

5. StoU, S. and Schweiger, A. (2006) J. Magn. Reson. 178(1), 42-55. 

6. Hill, S., Murugesu, M., and Christou, G. (2009) Phys. Rev. B 80, 174416. 

7. Quddusi, H. M., Liu, J., Singh, S., Heroux, K. J., delBarco, E., Hill, S., and 
Hendrickson, D. N. (2011) Phys. Rev. Lett. 106, 227201. 

8. Liu, J., Beedle, C. C., Quddusi, H. M., delBarco, E., Hendrickson, D. N., and 
Hill, S. (2011) Polyhedron 30, 2965-2968. 

9. Heroux, K. J., Quddusi, H. M., Liu, J., O'Brien, J. R., Nakano, M., delBarco, 
E., Hill, S., and Hendrickson, D. N. (2011) Inorganic Chemistry 50, 7367 7369. 

10. Beedle, C. C., Wang, W.-G., Koo, C., Liu, J., Zhou, A.-J., Nakano, M., OBricn, 
J. R., Wernsdorfer, W., Hill, S., Tong, M.-L., Chen, X.-M., and Hendrickson, 
D. N. (2013) Inorganic Chemistry X, XX-XX Under review. 

II. Hill, S., Pcrcnboom, J. A. A. J., Dalai, N. S., Hathaway, T., Stalcup, T., and 
Brooks, ,L S. (1998) Phys. Rev. Lett. 80, 2453-2456. 

12. Barra, A. L., Gatteschi, D., and Sessoli, R. (2000) Chemistry - A European 
Journal 6{9), 1608-1614. 

13. Takahashi, S., Edwards, R. S., North, J. M., Hill, S., and Dalai, N. S. (2004) 
Physical Review B 70(9), 094429. 

14. Kirman, C, Lawrence, J., HiU, S., Yang, E.-C, and Hendrickson, D. N. (2005) 
Journal of Applied Physics 97(10), 10M501. 

15. Carretta, S., Liviotti, E., Magnani, N., Santini, P., and Amoretti, G. (2004) 
Physical Review Letters 92(20), 207205. 

16. Wilson, A., Lawrence, J., Yang, E. C, Nakano, M., Hendrickson, D. N., and 
Hill, S. (2006) Physical Review B 74(14), 140403. 

17. Ramsey, C. M., delBarco, E., HiU, S., Shah, S. J., Beedle, C. C, and Hendrick- 
son, D. N. (2008) Nature Physics 4, 277-281. 

18. Hill, S., Datta, S., Liu, J., Inglis, R., Milios, C. J., Feng, P. L., Henderson, J. J., 
delBarco, E., Brechin, E. K., and Hendrickson, D. N. (2010) Dalton Transactions 
39(20), 4693-4707. 

19. Yang, E.-C, Kirman, C, Lawrence, J., Zakharov, L. N., Rheingold, A. L., Hill, 
S., and Hendrickson, D. N. (2005) Inorganic Chemistry 44(11) , 3827-3836. 



1 A Microscopic and Spectroscopic View of QTM 



27 



20. Maurice, R., deGraaf, C, and Guihery, N. (2010) Phys. Rev. B 81, 214427. 

21. Henderson, J. J., Koo, C, Feng, P. L., delBarco, E., Hill, S., Tupitsyn, I. S., 
Stamp, P. C. E., and Hendrickson, D. N. (2009) Physical Review Letters 103, 
017202. 

22. Feng, P. L., Koo, C, Henderson, J. J., Nakano, M., Hill, S., delBarco, E., and 
Hendrickson, D. N. (2008) Inorganic Chemistry 47 (19) , 8610-8612. 

23. Feng, P. L., Koo, C, Henderson, J. J., Manning, P., Nakano, M., delBarco, E., 
Hill, S., and Hendrickson, D. N. (2009) Inorganic Chemistry 48(8), 3480-3492. 

24. Redler, G., Lampropoulos, C., Datta, S., Koo, C, Stainatatos, T. C., Chakov, 
N. E., Christou, G., and Hill, S. (2009) Physical Review B 80(9), 094408. 

25. Koster, G. F., Dimmock, J. O., Wheeler, P. G., and Statz, H. (1963) Properties 
of the Thirty-Two Point Groups, M.l.T. Press, Cambridge, MA. 

26. Cotton, F. A. (1990) Chemical Apphcations of Group Theory, A Wiley- 
Interscience Publication, third edition. 

27. Liviotti, E., Carretta, S., and Amoretti, G. (2002) The Journal of Chemical 
Physics 117(7), 3361-3368. 

28. Barra, A.-L., Caneschi, A., Cornia, A., Gatteschi, D., Gorini, L., Heiniger, L.- 
P., Sessoli, R., and Sorace, L. (2007) Journal of the American Chemical Society 
129(35), 10754 10762. 

29. vanSlageren, J., Vongtragool, S., Gorshunov, B., Mukhin, A., and Dressel, M. 
(2009) Physical Review B 79, 224406. 

30. Wernsdorfer, W., Bhaduri, S., Boskovic, C., Christou, G., and Hendrickson, 

D. N. (2002) Physical Review B 65, 180403. 

31. Maccagnano, S., Acliey, R., Negusse, E., Lussier, A., Mola, M. M., Hill, S., and 
Dalai, N. S. (2001) Polyhedron 20, 1441. 

32. Park, K., Novotny, M. A., Dalai, N., Hill, S., and Rikvold, P. (2001) Physical 
Review B 65, 014426. 

33. Hill, S., Maccagnano, S., Park, K., Achey, R. M., North, J. M., and Dalai, N. S. 
(2002) Phys. Rev. B 65, 224410. 

34. Park, K., Novotny, M. A., Dalai, N. S., Hill, S., and Rikvold, P. (2002) Physical 
Review B 66, 144409. 

35. Chudnovsky, E. M. and Garanin, D. A. (2001) Physical Review Letters 87(18), 
187203. 

36. Chudnovsky E. M. and Garanin, D. A. (2002) Phys. Rev. B 65, 094423. 

37. Mertes, K. M., Suzuki, Y., Sarachik, M. P., Paltiel, Y., Shtrikman, H., Zeldov, 

E. , Rumberger, E. M., Hendrickson, D. N., and Christou, G. (2001) Physical 
Review Letters 87(22), 227205. 

38. delBarco, E., Kent, A. D., Rumberger, E. M., Hendrickson, D. N., and Christou, 
G. (2002) EuroPhys. Lett. 60, 768-774. 

39. delBarco, E., Kent, A. D., Rumberger, E. M., Hendrickson, D. N., and Christou, 
G. (2003) Phys. Rev. Lett. 91, 047203. 

40. delBarco, E., Kent, A. D., Hill, S., North, ,1. M., Dalai, N. S., Rumberger, E. M., 
Hendrickson, D. N., Chakov, N., and Christou, G. (2005) J. Low. Temp. Phys. 
140, 119-174. 

41. Cornia, A., Sessoli, R., Sorace, L., Gatteschi, D., Barra, A. L., and Daiguebonne, 
C. (2002) Phys. Rev. Lett. 89, 257201. 

42. Hill, S., Edwards, R. S., Jones, S. I., North, J. M., and Dalai, N. S. (2003) Phys. 
Rev. Lett. 90, 217204. 

43. Bircher, R., Cliaboussant, G., Sieber, A., Giidel, H. U., and Mutka, H. (2004) 
Phys. Rev. B 70, 212413. 



28 



Junjie Liu, Enrique del Barco, and Stephen Hill 



44. Hill, S., Anderson, N., Wilson, A., Takahashi, S., Petukhov, K., Chakov, N. E., 
Murugesu, M., North, J. M., dolBarco, E., Kent, A. D., Dalai, N. S., and Chris- 
tou, G. (2005) Polyhedron 24, 2284-2292. 

45. Chakov, N. E., Lee, S.-C, Harter, A. G., Kuhns, P. L., Reyes, A. P., Hill, S. O., 
Dalai, N. S., Wernsdorfer, W., Abboud, K. A., and Christou, G. (2006) J. Am. 
Chem. Soc. 128, 6975-6989. 

46. Subedi, P., Kent, A. D., Wen, B., Sarachik, M. P., Ycshurun, Y., Millis, A. J., 
Mukherjee, S., and Christou, G. (2012) Phys. Rev. B 85, 134441. 

47. Lampropoulos, C, Murugesu, M., Harter, A. G., Wernsdofer, W., Hill, S., Dalai, 
N. S., Abboud, K. A., and Christou, G. (2013) Inorg. Chem. 52, 258-272. 

48. delBarco, E., Kent, A. D., Chakov, N. E., Zakharov, L. N., Rhcingold, A. L., 
Hendrickson, D. N., and Christou, G. (2004) Phys. Rev. B 69, 020411. 

49. Carbonera, C, Luis, F., Campo, J., Sanchez-Marcos, J., Camon, A., Chaboy, 
J., Ruiz-Molina, D., Imaz, L, vanSlageren, J., Dengler, S., and Gonzalez, M. 
(2010) Phys. Rev. B 81, 014427. 

50. Burzurf, E., Carbonera, C, Luis, F., Ruiz-Molina, D., Lampropoulos, C, and 
Christou, G. (2009) Phys. Rev. B 80, 224428. 

51. Fernandez, J. F., Luis, F., and Bartolome, J. (1998) Phys. Rev. Lett. 80, 5659- 
5662. 

52. Luis, F., Bartolome, J., and Fernandez, J. F. (1998) Phys. Rev. B 57, 505-513. 

53. Prokof'ev, N. V. and Stamp, P. C. E. (1998) Phys. Rev. Lett. 80, 5794-5797. 

54. Lawrence, ,J., Yang, E.-C, Edwards, R., Olmstead, M. M., Ramsey, C, Dalai, 
N. S., Gantzel, P. K., Hill, S., and Hendrickson, D. N. (2008) Inorg. Chem. 47, 
1965-1974. 

55. Garg, A. (1993) EPL (Europhysics Letters) 22(3), 205. 

56. Park, C.-S. and Garg, A. (2002) Physical Review B 65, 064411. 

57. Li, F. and Garg, A. (2011) Physical Review B 83(13), 132401. 

58. Liu, J., delBarco, E., and Hill, S. (2012) Physical Review B 85(1), 012406. 

59. Yang, E.-C, Wernsdorfer, W., Hill, S., Edwards, R. S., Nakano, M., 
Maccagnano, S., Zakharov, L. N., Rheingold, A. L., Christou, G., and Hen- 
drickson, D. N. (2003) Polyhedron 22, 1727-1733. 

60. Yang, E.-C, Wernsdorfer, W., Zakharov, L. N., Karaki, Y., Yamaguchi, A., 
Isidro, R. M., Lu, G.-D., Wilson, S. A., Rheingold, A. L., Ishimoto, H., and 
Hendrickson, D. N. (2006) Inorg. Chem. 45, 529-546. 

61. Lawrence, J., Hill, S., Yang, E.-C, and Hendrickson, D. N. (2009) Phys. Chem. 
Chem. Phys. 2009(11), 6743-6749. 

62. Inglis, R., Taylor, S. M., Jones, L. F., Papaefstathiou, G. S., Perlepes, S. P., 
Datta, S., Hill, S., Wernsdorfer, W., and Brechin, E. K. (2009) Dalton Transac- 
tions 2009, 9157-9168. 

63. Inglis, R., Jones, L. F., Milios, C. J., Datta, S., Collins, A., Parsons, S., Werns- 
dorfer, W., Hill, S., Perlepes, S. P., Piligkos, S., and Brechin, E. K. (2009) Dalton 
Transactions 2009, 3403-3412. 

64. delBarco, E., Hill, S., Becdlc, C. C, Hendrickson, D. N., Tupitsyn, L S., and 
Stamp, P. C. E. (2010) Phys. Rev. B 82, 104426. 

65. Wernsdorfer, W., Stamatatos, T. C, and Christou, G. (2008) Phys. Rev. Lett. 
101, 237204. 

66. delBarco, E., Hill, S., and Hendrickson, D. N. (2009) Phys. Rev. Lett. 103, 
059701. 



1 A Microscopic and Spectroscopic View of QTM 



29 



67. Milios, C. J., Vinslava, A., Wernsdorfer, W., Moggach, S., Parsons, S., Perlepes, 
S. P., Christou, G., and Brechin, E. K. (2007) Journal of the American Chemical 
Society 129(10), 2754-2755. 

68. Waldmann, O. (2007) Inorg. Chem. 46, 10035 10037. 

69. Harman, W. H., Harris, T. D., Frcedman, D. E., Fong, H., Chang, A., Rinehart, 
J. D., Ozarowski, A., Sougrati, M. T., Grandjean, F., Long, G. J., Long, J. R., 
and Chang, C. J. (2010) J. Am. Chem. Sac. 132, 18115-18126. 

70. Zadrozny, J. M., Liu, J., Piro, N. A., Chang, C. J., Hill, S., and Long, J. R. 
(2012) Chem. Commun. 48, 3927-3929. 

71. Ruamps, R., Maurice, R., Boggio-Pasqual, M., Guihery, N., Batchelor, L., Liu, 
J., Hill, S., Mallah, T., and Barra, A.-L. (2001) J. Am. Chem. Soc. X, XX-XX 
doi:10.1039/B925358K (Article ASAP). 

72. Ishikawa, N., Sugita, M. Ishikawa, T., Koshihara, S., and Kaizu, Y. (2003) J. 
Am. Chem. Soc. 125, 8694-8695. 

73. AlDamen, M., Clomcntc-Juan, J. M., Coronado, E., Marti-Gastaldo, C, and 
Gaita-Ariho, A. (2008) J. Am. Chem. Soc. 130, 8874-8875. 

74. Rinehart, J. D. and Long, J. R. (2009) Journal of the American Chemical Society 
131, 12558-12559. 

75. Ghosh, S., Datta, S., Friend, L., Cardona-Serra, S., Coronado, E., and Hill, S. 
(2012) Dalton Trans. 41, 13697. 



Table 1.1. Comparison of tunneling gaps obtained for Mns from the MS and GSA 
models for resonances A; = 0, 1, 2 and 3, for the two cases ^ = (top) and 6 = 8.5° 
(bottom). 

fc n A GSA-gap (K) MS-gap (K) Ratio 

JT-axes parallel to the molecular z-axis 
2 Ag g 2.60 X 10-2 2.66 x lO'^ 0.98 
2 Ag^g 1-10 X IQ-^ 1.05 X lO-'^ 1.05 
2 2 A2,4 2.37 X 10-2 2.35 x 10"^ 1.01 



JT-axes tilted ^ = 8.5° away from the molecular z-axis 
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Fig. 1.1. Potential energy surface corresponding to the 2"'^-order anisotropy tensor. 
The surface is generated employing Eqn. 1.2 with \E/D\ — 1/5 and D < 0. The 
radial distance to the surface represents the energy of a spin as a function of its 
orientation. 




Fig. 1.2. Potential energy surfaces corresponding to the fourth order Stevens op- 
erators O4 (a), O4 (b), O4 (c) and O4 (d). As can be seen, the g-even operators 
include the xy-p\ane as an extra symmetry element while the g-odd operators in- 
clude improper rotations (see main text). 




Fig. 1.3. The molecular structure (a) and schematic representation of the magnetic 
core (b) of the Mng SMM. Color code: Mn — purple, Zn = green, O = red, N = 
blue, C — black and CI — dark gold. H-atoms have been omitted for clarity. 




Fig. 1.4. Zeeman diagram for a spin S — 6 multiplet with easy-axis anisotropy 
(73 < in Eqn. 1.8) and H/ /z. All possible non-zero tunneling gaps for C3 symmetry 
are labeled according to the scheme discussed in the main text. The inset shows the 
Ht dependence of the odd-n tunneling gaps. 



BlOl (Q„ symmetiy) 
JT-axes of Mn'" ions // the molecular C3 axis 




B\d\ (S^ symmetiy) 
JT-axes of Mn'^' ions tilted away from the molecular C3 axis 



Fig. 1.5. The influence of the orientations of the JT-axes of Mn ions on the zero- 
field magneto symmetry of the Mn^ SMM. In (a), the JT-axes of the Mn''^ ions 
(left) are parallel to the molecule C3 axis; consequently, the resultant Hamiltonian 
of the molecule (right) possesses Ceh symmetry. In (b), the JT-axes of the Mn^^' 
ions (left) are tilted away from the molecular C3 axis; consequently, the resultant 
Hamiltonian of the molecule possesses 5*6 symmetry 
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Fig. 1.6. Calculated QTM gaps for the Mna SMM as a function of the coupling 
constant J. Simulations were performed with the JT-axes tilted 8.5° away from the 
molecular Ca-axis. The QTM gaps associated with same |Am| value are rendered in 
the same color. Note that the results are plotted on a logarithmic scale. 




Fig. 1.7. (a) Field derivative of tlie magnetization curves obtained for a Mna single 
crystal at different temperatures, with B || z. (b) Ground-state tunnel splittings 
associated with resonances fc = (black squares), k — 1 (red circles), k — 2 (green 
triangles), and fc = 3 (blue stars) as a function of the transverse field Ht, with 
the JT-axes aligned along the C3 axis (thin lines) and tilted by 8.5 degrees away 
from the C3 axis (thick lines). The strength of the dipolar magnetic field in the 
sample is represented by the central gray area, with the corresponding splitting 
values achieved for such dipolar field values for resonances k = 1 and k = 2 (dashed 
horizontal lines). 
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Fig. 1.8. (a) The calculated BPI patterns associated with the ground state fc = 
QTM resonance for the Mua SMM: the color contour plot shows Ag g as function 
of Ht (with _B| set to zero); a compensating //l field is required that alternates 
between positive (red) and negative (blue) values, (b) The compensating //l field 
for Ag 5, as a function of the magnitude of Ht; note the curvature (except for the 
30° trace, for which Hi, = 0). 




Fig. 1.9. The structure (a) and schematic representation of the magnetic core (b) 
of the Ni4 SMM. Color code: Ni = oHve, O = red, N = blue, C = black and CI = 
dark gold. H-atoms have been omitted for clarity, (c) Representation of the zero- 
field magneto symmetry of the Ni4 SMM resulting from the situation in which the 
2"'^-order single-ion zfs tensors have their C2 axes tilted away from the molecular 
Si axis. Once added, the time reversal symmetry of the SO interaction guarantees 
that the resultant zero-field Hamiltonian of the molecule possesses Cih symmetry 
(see text for details). 
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Fig. 1.10. (a) Frequency dependence of the positions (in field) of the three EPR 
transitions associated with the isolated Ni'' triplet {S = 1) state in a diluted single 
crystal of the compound [Zn3.giNio.og(hmp)4(dmb)4Cl4] (see Ref. [19] for assign- 
ments of the A, B and C peaks). The colored curves correspond to best fits to 
the data employing the following single-ion zfs parameters: d — —5.30(5) cm~^, 
e = ±1.20(2) cm~^, — 2.30(5), and a tilting of the local z-axes of 15° away from 
the symmetry (c-) axis of the crystal. The energy splittings labeled A provide a 

direct measure of the tilting of the local d tensors; the dashed curves correspond 
to the non-tilted case, for which these splittings are zero. The widths of the col- 
ored curves reflect the uncertainty in the orientations of the local x- and y-axes, 
which were subsequently deduced from two axis rotation studies [19]. (b) 2D color 
map of the EPR absorption intensity as a fucntion of the magnetic field strength 
and its orientation within the hard plane of a single crystals of the S — 4 SMM 
[Ni4(hmp)(dmb)Cl]4 (see Ref. [61] for explanation of the peak labeling). Superim- 
posed on the absorption maxima (darker red regions) are fits (white curves) to the 
data that involve just a single adjustable zfs parameter, Bf = 4 x 10~* cm~^ (over 
and above those deduced from easy-axis measurements [16]). 
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Fig. 1.11. Zeeman diagram for the ground state 5 = 4 multiplet associated with 
the Ni4 SMM, simulated employing Eqn. 1.9. The fc = to 4 ground state QTM 
splittings are labeled in the figure. 




Fig. 1.12. The ground state QTM gaps for the Ni4 SMM as a function of Ht- 
The simulations were performed employing Eqn. 1.9 with the parameters given in 
the main text. The solid lines were generated with 6 = and the dash lines were 
generated with 9 = 15°. 




Fig. 1.13. The effect of disorder on the ground state QTM gaps for the Ni4 SMM. 
The simulations were performed employing Eqn. 1.9 by misaligning the zfs tensor of 
one of the Ni'^ ions with respect to the unperturbed molecular z-axis (see details in 
the main text). 




Fig. 1.14. (a) The Mn4-anca molecule, (b) The Mn4-Bet molecule. Color code: Mn 
= purple, O = red, N = blue, C = grey, H = white, B = pink and CI = green, (c) 
Sketch showing the different exchange interactions used to solve the four spin MS 
Hamiltonian for these molecules, (d) Trimer model representing the Mn4 molecules 
assuming an infinite J coupling between the two central Mn^'' ions. 
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Fig. 1.15. EPR spectra obtained at different temperatures with the field along the 
easy-axis at: 165 GHz for Mn4-anca (a); and 139.5 GHz for Mn4-Bet (b). Zeeman 
diagrams depicting the low- lying energy levels for Mn4-anca (c) and Mn4-Bet (d), 
obtained by diagonalization of the MS Hamiltonian (Eqn 1.6) using the trimer model 
depicted in Fig. 1.14d. Arrows relate the QTM peaks observed in the field deriva- 
tives of the magnetization curves obtained at different temperatures (bottom of the 
graphics) with the corresponding crossings between spin levels; with black arrows 
indicating crossings of the ground state 15 = 9, — 9) (red for excited states) with 
other states within the same spin multiplet {S — 9), and blue arrows signaling both 
ground and excited crossings involving levels of different spin length (5' = 9 and 
5 = 8). 
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Fig. 1.16. Modulation of the QTM probabilities for resonances fc = (a) and 
k = 1 (b) as a function of Ht applied at different angles, cj>, within the xy-plane of 
the Mn4-Bet SMM. The asymmetry of the BPI pattern of oscillations in resonance 
A: = 1 is inverted upon reversal of Hl- (c) and (d) show the modulation of the tunnel 
splitings of resonances = and k = 1, respectively, for different directions of the 
transverse field. 



